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Preface 



The conjúgate nature is an inherent feature of any heat transfer problem, 
because the intention of any heat transfer process is an interaction of at least 
two médiums or subjects. This inherent property specifies the fact that heat 
transfer processes began to be studied as conjúgate problems soon after such 
a possibility was realized with the advent of computers. Apparently, the 
same consideration gives an understanding why the number of publications 
in this area continúes to grow like an avalanche. 

The development of the conjúgate heat transfer approach is important not 
only in that conjúgate formulation is one of fundamental nature access, but 
also that it has a wide range of applications. The advanced conjúgate model- 
ing of convective heat transfer phenomena is now used extensively in differ- 
ent applications. Beginning with simple examples in 1965 through the 1970s, 
this approach is now used to create models of various device operational 
and technological processes from relatively simple procedures to complex, 
multistage, nonlinear processes. 

I have studied the conjúgate convective heat transfer systematically for 
more than 40 years, beginning in about 1965, and initially this book was 
intended as a sum of my achievements. However, in the process of work- 
ing, it became clear that a real valué of any result can be understood only 
from the point of view of the contemporary knowledge in the area of inter- 
est, and the book gradually has been transformed into a review of the body 
of publications. Thus, finally, the book has become one of the first consistent 
presentations of the current situation in conjúgate convective heat transfer. I 
hope that I have been sufficiently successful in meeting and resolving such 
a challenge. 

This book begins with a short introduction, "What is conjúgate heat trans- 
fer like?" and closes by discussing the question, "Should any heat transfer 
problem be considered as a conjúgate?" Three parts and ten chapters incor- 
pórate the theory, the methods of study, and the applications of conjúgate 
convective heat transfer. More than a hundred examples of solved typical 
conjúgate problems in different areas from early publications to many con- 
temporary results are considered. Although it is obvious that the choice of 
such examples is random and depends on preferences, background, and 
insight, I hope that the presented pattern of examples offers the reader a 
basic understanding of the current situation in the area of conjúgate convec- 
tive heat transfer. 

The first part of the book consists of the approximate analytical Solutions 
of conjúgate heat transfer problems (Chapter 2) and methods of estimation of 
heat transfer from nonisothermal surfaces (Chapter 1), which were basically 
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used in the early works to get approximate conjúgate Solutions. Two types of 
approximate Solutions are considered. The first ineludes Solutions based on 
the approximate energy equation when the approximate velocity distribu- 
tion across the thermal boundary layer is used. Another group contains Solu- 
tions with results that could not be refined or that are difficult to upgrade. 
Most examples of this type of solution are taken from early works, yet some 
results obtained recently are included as well, as these approximate methods 
are still in use due to their simplicity and comprehensibility. 

This first part of the book, presents the initial period of conjúgate heat 
transfer development, and the remainder ineludes the theory, the methods 
of Solutions, the applications, and the examples that fit the modern level of 
knowledge in this area. 

The second part contains five chapters (Chapters 3 through 7). These chap- 
ters represent theory (Chapters 3 through 5); modern analytical (Chapter 6); 
and numerical (Chapter 7) methods of conjúgate convective heat transfer 
investigation and corresponding examples of Solutions. Chapters 3 and 4 
incorpórate the exact and accurate approximate Solutions of the thermal 
boundary layer equations for arbitrary nonisothermal surfaces in the lami- 
nar (Chapter 3) and turbulent (Chapter 4) flows. Solutions for both the lami- 
nar and turbulent incompressible flows with zero and nonzero pressure 
gradients are given in two forms: in the differential form in the series of 
consecutive derivatives of the temperature; and in the integral form with 
the influence function of the unheated zone. In the case of laminar flow, the 
other exact Solutions for arbitrary nonisothermal surfaces are obtained: for 
the cases of heat flux distribution and unsteady temperature distribution; 
for high-speed compressible flow; for power law non-Newtonian fluids; for 
fluids with significant energy dissipation; and for axisymmetric and rotation 
bodies. 

The general theory of nonisothermal and conjúgate heat transfer devel- 
oped on the basis of the Solutions presented in Chapters 3 and 4 is outlined 
in Chapter 5. The general properties of nonisothermal and conjúgate convec- 
tive heat transfer are formulated, and the basic parameters determining the 
intensity of heat transfer are indicated. The effect of different factors is stud- 
ied, including the temperature head distribution, the temperature head and 
pressure gradients, the flow regime, the body shape, and the Biot number. 
Reynolds analogy and gradient analogy (between pressure and temperature 
head gradients) are investigated, and conditions are formulated when such 
analogies are valid. A phenomenon of heat flux inversión that is similar to 
separation of flow is studied, and the like and unlike properties of both phe- 
nomena are indicated. The deformation of the temperature head profile in 
the thermal boundary layer as a physical reason for the existence of heat 
inversión is analyzed. On the basis of this analysis, it is shown that there are 
temperature head distributions that correspond to surfaces with zero heat 
transfer. Such a surface can be theoretically generated, if the temperature 
head profile has a zero angle tangent at each point of the surface. In the last 
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section of this chapter, the technique of optimization of heat transfer in flows 
over bodies is provided. 

In Chapter 6, the general method for solving conjúgate heat transfer prob- 
lems is developed. Both the differential and integral forms for surface heat 
flux obtained in Chapter 3 are used. It is shown that the differential expres- 
sion for heat flux is, in fact, the general condition of the third kind, which 
transforms in the usual condition of the third kind for the case of an isother- 
mal surface when all of the derivatives become zero. By taking into account 
some derivatives, a more accurate solution can be obtained, and it becomes 
possible to estimate an error, which can occur when the condition of the 
third kind is used, and then to decide whether there is a need for a conjúgate 
problem solution. The method of reducing the steady or unsteady conjúgate 
heat transfer problems to an equivalent problem for a conduction equation 
with a general boundary condition of the third kind is offered, which is simi- 
lar to a well-known approach with the usual condition of the third kind. 
Using both the differential and integral expressions for surface heat flux, 
the method of successive differential-integral approximations is suggested 
in which integral expression is used to refine the result obtained by the dif- 
ferential expression, and vice versa. The simple method for calculating the 
conjúgate heat transfer from thermally thin bodies using universal tabulated 
functions is proposed for both the steady and unsteady cases. It is shown that 
the singularity in temperature distribution at the leading edge of the píate 
is defined by the Reynolds number exponent in the Nussel number equa- 
tion for an isothermal surface, so that the series describing the temperature 
distribution at the leading edge is in integer powers of one-half and one-fifth 
for laminar and turbulent flows, respectively. The other cases are considered, 
and a general rule for determining the character of singularity in the tem- 
perature distribution at the leading edge is derived. It is also proved that the 
Biot number is a conjúgate criterion defining the level of thermal coupling 
between the fluid and the solid, and that the various relations suggested by 
different authors are in fact the other forms of the Biot number. 

Sections 6.7 through 6.10 in Chapter 6 discuss the following analytical 
methods currently used for solving conjúgate heat transfer problems: inte- 
gral transforms and similar methods, asymptotic series in eigenfunctions, 
superposition, Green's function, and perturbation methods. Examples are 
presented for each method. 

Chapter 7 contains an analysis of contemporary numerical methods for 
solving conjúgate heat transfer problems. The relation and mutual impor- 
tance of numerical and analytical approaches are discussed. Different finite- 
difference and finite-element methods are described and compared using the 
weighted residual approach. The difficulties in computing the convection- 
diffusion terms, velocities, and pressure, and methods of resolving these 
problems are discussed. The central-difference and upwind scheme as well 
as some others are compared to show their advantages and shortcomings. 
The SIMPLE and SIMPLER software are described. The existing numerical 
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approaches for subdomain Solutions conjugation, such as simultaneous solv- 
ing of one large set of equations for fluids and solids, iterative schemes, and 
using the superposition principie for conjugation are analyzed. Some mod- 
ern numerical approaches to improve the convergence of the iterative conju- 
gation methods, such as the multigrid or meshless techniques, are discussed 
as well. A number of examples of the numerical Solutions of the conjúgate 
heat transfer of the flows in tubes, channels, and around and inside the bod- 
ies are presented to illustrate the use of the described methods. 

The third part of the book is composed of three chapters describing appli- 
cations in the thermal treatment of material (Chapter 8), in the technologi- 
cal processes (Chapter 9), and in the manufacturing equipment operation 
(Chapter 10). Numerous conjúgate heat transfer problems of simulation of 
various industrial processes are presented including: 

• Protection of the reentry rocket by a thermal shield (Example 10.10). 

• Heat transfer from a moving píate (Examples 8.1 and 8.2). 

• The thermal process in the microchannel of challenging a heat 
exchanger (Example 10.7). 

• The Czochralski crystal growth process (Example 9.4). 

• The optical fiber coating process (Example 8.7). 

• The rewetting process in nuclear reactors in case of an emergency 
accident (Examples 10.11 and 10.12). 

• The cooled turbine blades (Examples 10.8 and 10.9). 

• The multidimensional freeze-drying process of food (Example 9.9). 

Different approaches are used. Mostly, the finite-difference and finite-element 
methods are employed because the problems are usually complicated. 
However, some analytical and other numerical methods are used as well. In 
particular, these inelude: 

• The special method of separation variables is used in the problem of 
rewetting the surfaces by falling film (Example 10.11). 

• The integral method is employed in the problem of rewetting semi- 
infinite piafes (Example 10.12). 

• The method of reduction of the conjúgate heat transfer problem to 
an equivalent heat conduction problem (Chapter 3) is applied in the 
problem of drying a moving strip of material (Example 8.2). 



Abram Dorfman 
Ann Arbor, Michigan 
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Introduction 



What Is Conjúgate Heat Transfer Like? 

Exact formulation of convective heat transfer problems requires using the 
boundary condition of the fourth kind. Such a boundary condition consists 
of the conjugation of the body and fluid temperature fields at their interface. 
Because a conjugation procedure is complicated, instead of it, the boundary 
condition of the third kind has been used since the time of Newton: 

q m = h (T w -TJ ( 1 . 1 ) 

The heat transfer coefficient, h, is usually determined experimentally, and 
no well-found theoretical approach was available until the last few decades. 
Because the valué of the heat transfer coefficient depends on the wall tem- 
perature distribution, data for the isothermal wall (T = const ) heat transfer 
coefficient have been used for practical calculations. This simplified approach 
in which the effect of actual wall temperature distribution is neglected was 
acceptable before the advent of the Computer, when the required calculation 
accuracy was not as high. 

Intense interest in the conjúgate convective heat transfer aróse at the end 
of 1960s, and beginning at that time, many convective heat transfer problems 
have been considered using a conjúgate, coupled, or adjoint formulation. 
These three equivalent terms correspond to a situation when the solution 
domain consists of two or more subdomains in which studied phenomena 
are described by different differential equations. After solving the problem 
in each of the subdomains, these Solutions should be conjugated. The same 
procedure is needed if the problem is governed by one differential equation, 
but the subdomains have different materials or other properties. 

For example, heat transfer between a body and a fluid flowing past it is 
a conjúgate problem, because the heat transfer inside the body is governed 
by the elliptic Laplace equation or by the parabolic differential equation, 
while the heat transfer inside the flowing fluid is governed by the elliptic 
Navier-Stokes equation or by the parabolic boundary layer equation. The 
solution of such a problem gives the temperature and heat flux distributions 
along the body-fluid interface, and there is no need for a heat transfer coef- 
ficient. Moreover, the heat transfer coefficient can be calculated using these 
results. Another example of a conjúgate problem is the transient heat transfer 
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between a hot píate and a cooling thin fluid film flowing along it. In such a 
case, the píate at each moment is divided into a wet part covered with mov- 
ing film and a dry one until uncovered. Although heat transfer in both por- 
tions of the píate is governed by the same conduction equation, the Solutions 
obtained for each of the parts should be conjugated, because the thermal 
properties of wet and dry parts are different. 

There are many conjúgate problems in other areas of Science. For instance, 
studying subsonic-supersonic flows required conjugation, because the sub- 
sonic flow is governed by elliptic or parabolic differential equations, while 
the supersonic flow is described by the hyperbolic differential equation [1]. 
Combustión theory and biological processes are two other examples. Every 
combustión process has two areas containing fresh and burnt gases with dif- 
ferent properties [2]. In biology, diffusion processes usually proceed simul- 
taneously in qualitatively different areas (e.g., in membranes) and, therefore, 
require a conjugation procedure [3]. 

Analogous mathematical problems with a similar formulation, usually 
called mixed problems, began to be considered much earlier than these 
physical Systems. The most famous mixed parabolic/hyperbolic equation 
was investigated by Tricomi in 1923. This equation is now widely used for 
subsonic-supersonic flows [1]. 
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Nomenclature 




í) A 

Br = — — Pr” ! Re" 

u 



c c 

'-'I' ^2 



c 



/ 




Biot number 



Brun number 

Exponents of influence functions in integral 
forms determining surface heat flux and 
temperature 

Friction coefficient 



2St/Cf 



Reynolds analogy coefficient 



c,c 



p 



Specific heat and specific heat at constant 
pressure J/kg K 



c = pcA 



Thermal capacity J/m 2 K 



D, D h 
D m 



Wx) 



Diameter and hydraulic diameter m 
Diffusion coefficient, m 2 /s 

Eckert number 

Influence function of unheated zone at tem- 
perature jump 



f,G/x) 



Influence function of unheated zone at heat 
flux jump 




Fourier number 



G = 



Mi 
KA T 



Griffith number, Example 8.6 
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Nomenclature 


Gr 

V 2 


Grashof number 


Sk'k 


Coefficients of series in differential forms 
determining surface temperature and heat 
flux 


$ 


Gravitational acceleration, m/s 2 


kK 


Heat and mass transfer coefficients, W/m 2 K 


k 


Specific heat ratio or turbulence energy in 
k-co model 




Constants in power rheology laws for non- 
Newtonian fluids 


Kn = — 
Dh 


Knudsen number, Example 10.7 


1 


Body length or mixing length or free path, m 


L 


Characteristic length, m 


D 

Le = — — 
a 


Lewis number, Example 8.8 


A h 

Ls = - 

PfrlUl A 


Leidenfrost number, Example 10.12 


Lu- P C 

Pzv^w 


Luikov number, Section 6.6.1 


M = -^- 


Mach number 


M 


Moisture contení, kg/kg 


n, s 


Exponents in power rheology law for non- 
Newtonian fluids 


m kL hL s+1 

Nu=— ,Nu = 

A K U s 

‘ i 


Nusselt number for Newtonian and non- 
Newtonian fluids 


P 


Pressure, Pa 




Nomenclature 
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Peclet number for Newtonian and non- 
Newtonian fluids 



Pr = 





Prandtl number for Newtonian and non- 
Newtonian fluids 



Heat flux, W/m 2 , and volumetric heat source, 
W/m 3 

r/s Exponent in heat transfer coefficient expres- 

sion for isothermal surface 



Ra=«^ 

va 

Re= ÜL 

v K t 



Rayleigh number 

Reynolds number for Newtonian and non- 
Newtonian fluids 



Se = 



D,„ 



Sh = 



pc p D m 



St = 



Sk = 



Ste = 



h 

pc p U 

4 crTfL 

c A T 



_ _p_ 



A 



f 

T 

u, v,w 



U 



Schmidt number, Example 9.1 

Sherwood number, Example 9.1 

Stanton number 

Starks number, Example 7.15 

Stephan number, Example 9.2 
Time, s 

Temperature, K 

Velocity components; u, v parts in integrating 
by parts 

Velocity on outer edge of boundary layer 
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Nomenclature 



u t = \¡ z u. ! P Friction velocity 

u + = u/ u T , y + = \ju x / v Variables in wall law 

x,\j,z Coordinates 



Greek Symbols 

a 

P 



h K 

h h 



Xf 



c f . 






8,* h 



K 

A 

A 



M 

v 

z 

e=T-T 



Thermal diffusivity m 2 /s 

Dimensionless pressure gradient in self- 
similar Solutions and turbulent equilibrium 
boundary layer or volumetric thermal expan- 
sión coefficient, 1/K 

Nonisothermicity and nonisobaricity coeffi- 
cients, Example 8.8 

Nonisotachicity coefficient, Section 5.2 

Boundary layer thicknesses, m 
Delta function and Kronecker delta 
Body or wall thickness, m 
Constant determining mixing length 
Thermal conductivity, W/mK 
Latent heat, J/kg or AJA 
Viscosity, kg/s m 
Kinematic viscosity, m 2 /s 
Unheated zone length, m 
Temperature excess 




Nomenclature 
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e = T -T 

W W OO 

P 

a 

T 

®,<p 

V 

(0 



Temperature head 
Dcnsity, kg/m 3 

Stefan-Boltzmann constant, W/m 2 K 4 
Shear stress, N/m 2 
Prandtl-Mises-Górtler variables 
Stream function, m 2 /s 
Specific dissipation rate in k-co model 



Some of these symbols are also used in different ways, as indicated in each 
case. 



Subscripts 

av 

nd 

as 

bl 

e 

i 

L 

m 

o 

s 



Average 

Adiabatic 

Asymptotic 

Bulk 

End or effective or entrance 
Initial; inside 
At x = L 

Mass average, or mean valué, or moisture 

Outside 

Solid 



sd 



Sound 
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Nomenclature 



t 

tb 

w 

$ 



oo 



* 



Constant heat flux 

Thermal 

Turbulent 

Fluid-solid interface 
After jump 
Far from solid 
Isothermal 



Superscipts 

+, - From both sides of interface 

+, + + Wall law and other variables for turbulent 

boundary layer 



Overscores 

o, o 



Dimensionless; or transformed or auxiliary 




Part I 

Approximate Solutions 





Analytical Methods for the Estimation of 
Heat Transferfrom Nonisothermal Walls 



In this chapter, the approximate analytical methods for calculating heat 
transfer from arbitrary nonisothermal walls are analyzed. Most of the early 
Solutions of conjúgate heat transfer problems are based on those approximate 
methods. Moreover, approximate analytical methods are frequently used at 
present due to their simplicity and easy physical interpretation of the results. 



1.1 Basic Equations 

The majority of practically important problems of flow and heat transfer are 
characterized by high Reynolds (Re) and Peclet (Pe) numbers. In such a case, 
viscosity and conductivity are significant only in thin boundary layers, and 
the System of Navier-Stokes and energy equations is simplified to bound- 
ary layer equations. For laminar steady-state flow of an incompressible fluid 
with constant thermophysical properties, these equations are [1] 



du dv 
dx dy 



du du , T dU d 2 u 

U — h 17 — U V- 

dx dy dx dy 2 



= 0 , 



dT dT d 2 T 
u^- + v- — a 



r duX 



dx dy dy 2 c p {dy , 



= 0 



( 1 . 1 ) 



The Solutions of System (1.1) must satisfy the boundary conditions at the 
body surface and far away from the body: 

y = 0 u = v = 0, T = TJx), y->°° u -> U, T-^T, (1.2) 



Prandtl and Misis independently transformed the boundary layer equa- 
tions into the form that is cióse to that of the heat-conduction equation [1]. This 
form is obtained using new variables x, stream function y/, and Z = U 2 - u 2 : 



dy/ dy/ dZ d 2 Z 
U ~ldy' V ~~lh' dx~ VU Yy T 2 

y/ = 0 Z = U 2 (x), T = T(x), 



dT 
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f dT ) 


V j 


í du ) 
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u \ 




dx 


dy/ 


l ¿Y) 


1 c p 1 





y/ — / °° Z — / 0, T — / 
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Equation (1.3) has a singularity on the surface such that the derivatives with 
respect to y/ become infinite. Because the velocity near the wall is an analytic 
function, one gets 



u ~ c r y + c 2 y 2 + ., , 




c c 2 
~- 1 v 2 +— y 3 + - 
2 J 3 J 



(1.4) 



Henee, y ~ y/ 1/2 , the derivative dy / dy/ ~ y/~ l/2 , and dy/dy/ — > °° when V > 0- 



1.2 Self-Similar Solutions of the Boundary Layer Equations 

If the body has a shape of a wedge with opening angle np, then the veloc- 
ity of the potential of an ideal fluid flow obeys a power law distribu- 
tion. If the temperature head along the wedge also obeys a power law 
distribution, 

U(x)=Cx m , m = p/(2-p) r T w (x)-T x =C 1 x mi (1.5) 



and then the partial differential equation (Equation 1.1) reduces to ordinary 
differential equations. In this case, it is said that the boundary layer equa- 
tions have self-similar Solutions. Physically, this means that the velocity and 
the temperature head distributions in different sections of the boundary 
layer are identical, so that single velocity and single temperature head pro- 
files form if they are plotted in the self-similar variables: 



u 

U 






T-T 

w 

T -T 

°° W 



0(?7) r¡ = y 



(m + l)Cx m ~ 1 
2v 



(p = ¥ 



m + 1 
2vCx m+1 



( 1 . 6 ) 



Equation (1.1) in the variables (p,9, and i] transform into ordinary differential 
equations: 

(p"' + yxp" + p(l-(p' 2 ) = 0 9" + Pr(p9' + l -Fv(p'(l-9) = 0 (1.7) 

m + 1 



Equation (1.7) is numerically integrated and tabulated [1], Using these 
results, the friction coefficient, the Nusselt number, and the heat flux are 
determined: 



C/ = 



= &>, Nu,=^xo ,.=c, 



m - 1 

l- 

2 



( 1 . 8 ) 
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TABLE 1.1 

Self-Similar Solutions 



m + 1 

Thermal J ^ ® ®) 



Boundary 





Condition 


P 


m 


^2(m + l))<p"(0) 


m 1 


Pr— >0 


Pr = l 


Pr — »<*> 


Píate 


T w = const. 


0 


0 


0.664 


0 


0.564Pr 05 


0.332Pr 1/3 


0.339Pr 1/3 


Píate 


q = const. 


0 


0 


0.664 


1/2 


0.885Pr 05 


0.453Pr 1/3 


0.458Pr 1/3 


Stagnation T u , = const. 


n/2 


1 


2.466 


0 








point 


q w = const. 


k/2 


1 


2.466 


0 


0.791Pr 05 


0.570Pr 1/3 


0.661Pr 1/3 



It is seen that in self-similar flows, the heat flux changes along the wall also 
according to the power law with exponent m 2 = m 1 + (1 - m)/2. Thus, heat flux 
is constant on the wall if m, = (1 - m)/ 2. Henee, on the píate (m = 0) with con- 
stant heat flux, the temperature varies as ~x 1/2 , while the stagnation point 
surface (m = 1) in such a case is isothermal. 

Example 1.1: The Effect of Boundary Conditions 

Table 1.1 shows that boundary conditions significantly affect the valué 
of the Nusselt number. Thus, on the stagnation point surface, the heat 
transfer coefficients for the surfaces with T w = const. and q w = const. are 
equal, while for the píate these coefficients differ remarkably. For the 
médium and high Pr, the difference between two coefficients is 36%, 
whereas for the low Prandtl numbers it reaches 57%. 



1.3 Solutions of the Boundary Layer 
Equations in the Power Series 

Self-similar Solutions are applicable to only a narrow class of problems. In 
general, when the velocity U(x ) and temperature head T w - T r = 6fx) are arbi- 
trary f unctions, the boundary layer equations may be solved by approximate 
methods. In 1908, Blasius was the first to suggest such a method [1]. The 
method is based on the assumption that the velocity and temperature head 
distributions around the symmetrical body can be presented in the form of 
a power series: 

U(x) = u x x + u 3 x 3 + u 5 x 5 + 0 w (x) = 0. v(j + 0 w2 x 2 +0 w 4 x 4 + - (1-9) 

The Solutions presented in the same form of power law series. 



y/(x, y) = yq {y)x + y/ 3 (y)x 3 + iy 5 (ij)x 5 + ■■■ 
T-T m = d(x, y) = e 0 (y)+ 9 2 (y)x 2 + d 4 {y)x 4 + - 



( 1 . 10 ) 
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being substituted into boundary layer equations give the System of ordinary 
differential equations defining the coefficients of these series. Howarth [1] 
shows that the variable ri = yJu 1 /x transforms equations for series coeffi- 
cients into a universal form that is independent on a specific problem. These 
coefficients and corresponding derivatives are tabulated [2,3] and can be used 
to calcúlate the friction and heat transfer for any series (Equation 1.9). For air 
(Pr = 0.7), the Nusselt number is given as (D is a characteristic length) [3]: 



Nu = 



D I u 

íT 



w 0 



0.4959 + 0.4764 — x 2 + — 

u, IL 



1.054 + 0.6678- 



ll l U 3 J 



+ 0 ... 



0.852x 2 + 06678— x 4 + — 
u. 



( 1 . 11 ) 



+ 0 J1.054X 4 + •••] + • 



The more general case of an asymmetrical body when a power law series 
(Equation 1.9) contains all valúes of exponents is considered in the literature [3]. 

Computational experience shows that the methods based on the power 
law extensión yield satisfactory results basically at a relatively small distance 
from the stagnation point. In order to obtain satisfactory accuracy at a large 
distance from the stagnation point, a large number of universal functions 
is required. This is because the number of universal functions needed for 
expressing coefficients in final series increases rapidly with the number of 
terms. This problem is especially severe for the asymmetrical bodies. At the 
same time, in the simpler case of symmetrical flow over a body, universal 
functions for the shear stress have already been calculated up to terms con- 
taining x 11 , and even the velocity profiles on the circular cylinder including 
a point of separation are calculated [1]. 

Gortler [1] suggested power law series Solutions in special variables: 



0 = — 
v 



X 



yu 

0 — i — ' 
vy20 



(p = 



¥ 



'V20 ' 



o= 



T-T 

w_ 

T -T 



( 1 . 12 ) 



Converting Equation (1.1) to these variables, one obtains 



d 3 (p d 2 (p 



1- 



r d(p^ 



= 20 



d 2 (p di 'p d(p d 2 (p 



d<¡>dr¡ di) 30 drj 2 



(1.13) 
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di -p d9 
30 di] 
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U0 W dx 



(1.14) 
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If the velocity U(x) and the temperature head O m (x) = T w {x)-T r are power 
law functions (Equation 1.5), then /3 and /3, are constant. In this case, (p and 
9 are functions only of variable O, and the Gortler Equations (1.13) and (1.14) 
become identical to the self-similar Equation (1.7). In general, when U(x) and 
9 w (x) are arbitrary functions, /3 and /3, are presented as power series in vari- 
able O. In this case, the solution of Equations (1.13) and (1.14) can be introduced 
in the form of a similar power law series in O. Then, the ordinary differential 
equations for coefficients of the Gortler series are obtained, which are analo- 
gous to these for the Blasius series. The Gortler series are used for solving a 
number of problems. In particular, these series are used to study forced [3] 
and natural [4] convective heat transfer from nonisothermal surfaces. It is 
found that the Gortler series converge more rapidly than the Blasius power 
series, but their convergence is also not always satisfactory. 



1.4 Integral Methods 

Integral methods are widely used because of their relative simplicity, which 
follows from replacing the differential equations of the boundary layer by 
the integral relationships. In this case, the conservation laws are not satisfied 
in each point of the flow, but only integrally, on the average for each cross- 
section through the boundary layer. 

To obtain the integral relationships, Equation (1.1) is rewritten in the form 



3(1 lu) | á(lfo) dU 

dx 3 y dx 



3 (n 2 ) + 3 (nv) ^ dU 3 2 u 

dx 3 y dx 3 y 2 



(1.15) 



The first part of Equation (1.15) is obtained by multiplication of the continu- 
ity Equation (1.1) by U. The second part is a sum of continuity Equation (1.1), 
multiplied by u, and momentum Equation (1.1). Integrating across the bound- 
ary layer the difference between the first and second parts of Equation (1.15) 
and taking into account that the derivative (dn/dy) at y — > °° and the product 
v(U - u) at y = 0 and y — > °° turn to zero, one obtains the integral boundary 
layer momentum equation. The integral boundary layer energy equation is 
obtained in a similar manner by integrating across the boundary layer the 
difference between the energy Equation (1.1) and continuity Equation (1.1) 
multiplied by (T -TJ: 



— (U8 n ) + U— 8* = ^- 
dx dx p 



dx 



(U9JD 



P C P 



(1.16) 



8 * = 




n 

U 



dy 



8 ” = 




u 

U 



dy 



<?;* = 



f « t-t„ 



(1.17) 
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Here 8* ,8** ,8*" are displacement, momentum, and energy thicknesses. 
Because the integral boundary layer equations satisfied the conservation 
laws on the average for each cross-section through boundary layer, they 
do not yield any information on the velocity and temperature distributions 
across the boundary layer. For this reason, the velocity and temperature pro- 
files are found in integral methods by selecting suitable sets of functions 
satisfying the boundary conditions. 



Example 1.2: Friction and Heat Transfer from 
an Arbitrary Nonisothermal Píate [5] 

The velocity and temperature distributions are specified as the four 
power polynomials. Ten coefficients of these polynomials are selected 
in such a manner as to satisfy the following ten conditions on the píate, 
as well as on the outer edge of the velocity (y = 5) and thermal (y = 8 1 ) 



layers. 












y=o. 


u = 0, 


T = TJx); 


d 2 u d 2 T 

d y 2 " dy 2 ■ ' 






y =8, 


u = U, 


du d 2 u 
3y 9y 2 


II 

s» 

H 

II 


dT _ 
8y 


— = 0 
9y 2 



( 1 . 18 ) 



The third conditions at y = 0 follow from the boundary layer equation 
for gradientless flow, and the other conditions are obvious. Simple cal- 
culations give the profiles 

77 = 2 ? 7 - 2 tj 3 + r¡* = F(r]) r¡= y ^ = 2ri t - 2ijf + r¡f = r¡,= y 

u 8 T - T , ó, 



( 1 . 19 ) 



Then, Equation (1.16) yields the differential equations, defining the 
thicknesses of the velocity and thermal boundary layers. The dynamic 
thickness is obtained simply: 

= ( 1 . 20 ) 

Two cases must be considered to get the thermal boundary layer thick- 
nesses: when the thermal boundary layer is thinner (Pr > 1) and thicker 
Pr — > 0 than the velocity boundary layer. Taking into account that in the 
second case outside of the velocity boundary layer, the function F t (r\ t ) = 1, 
relations for the first and second cases are obtained (e = 8 t / 8): 



8 .** 2e 3e 3 e 4 , _ 5" 3 3 2 3 

— — 1 1 , (c ^ 1 ). — 1 

8 t 15 140 180 8 t 10 lOe 15e 2 140e 4 



1 

180e 5 



(£> 1 ) 



( 1 . 21 ) 
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Retaining here only the two first terms and putting the restricted 
expressions into the second Equation (1.20) gives two ordinary differ- 
ential equations. The solution of these equations yields two expressions 
defining the thickness of the thermal boundary layers for (Pr > 1) and 
Pr — t 0: 



8 t _ 0.871 

8 ~ Pr l/3 [d m (x)¡ ,/2 x í/4 






S t _ 0.626 

8 ~ Vx 1,z Q m (x)x in 




( 1 . 22 ) 



(1.23) 



Comparing these results with self-similar Solutions gives the estimation 
of accuracy. 

Nusselt numbers corresponding to Equation (1.8) for a power law tem- 
perature head are 



Nu r = 0.358(2/»! + 1) 1/3 Pr 1 / 3 Re 3/2 , Nu r = 0.548(27 ?/j + 1) 1 ' 2 Pr 1 / 2 Re 3/2 



(1.24) 



Then, for surfaces with 9 U1 = const. (/«j = 0) and q w = const. (///j = 1/2), one 
obtains: 0.358,0.548 and 0.452,0.775. for (Pr > 1) and Pr — > 0, respectively. 
According to data from Table 1.1, the largest difference, about 12%, is in 
the case of q w = const. and Pr — > 0. 

The case when both functions 6 ul (x ) and U(x) are arbitrary is examined 
in [6]. 

Integral methods employing polynomials generally yield satisfac- 
tory results in cases when the profiles of velocity and temperature in 
the boundary layer do not exhibit singularities and can be accurately 
described by polynomials. This condition generally holds at relatively 
small velocity gradients in the mainstream and small temperature heads. 
In the case of large, especially of large negative, gradients (divergent 
flows and a temperature head falling along the flow), approximation of 
velocity and temperature profiles in a boundary layer by polynomials 
becomes unsatisfactory. 

In these cases, significantly better results are obtained when the veloc- 
ity and temperature distributions across the layer are described by a set 
of profiles corresponding to the exact solution of some particular prob- 
lem of the boundary layer, for example, by a set of profiles obtained via 
self-similar Solutions. Such a method for the velocity boundary layer was 
developed by Kochin and Loytsyanskiy [7]. Subsequently, this method 
for heat transfer coefficients was developed by Skopets [8], Zysina- 
Molozhen et al. [9], and Drake [10]. 

The Kochin-Loytsyanskiy method is based on an assumption that 
at arbitrary U(x) and 0 w (x), the velocity and temperature distributions 
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across the layer are described by function (1.6) that depends on param- 
eters, just as in the case of self-similar flows: 

U z /u = <p\ n,p), (T - T w )/ (T m - T w ) = d(r¡, ¡i, Wj, Pr), r 1 = yp'/¡iv 

(1.25) 



The parameters ¡3 and m, are selected in such a manner that the inte- 
gral relation (1.16) is satisfied in each section across the boundary layer. 
For this purpose, function (1.25) is substituted into Equation (1.16). This 
yields two differential equations: 



dx U' 1 U J dx 



w e WJ 






(1.26) 



which define, at specified U(x) and 9 w (x), functions f(x) = 8 ,,2 U '/ v and 
f t (x) = S; ,2 U'/ v, termed the velocity and thermal shape parameters, 
respectively. Functions F(f) and F, (/,/,) are evaluated from tables of 
self-similar Solutions. For known F(f) and F f (/,/,), Equation (1.26) 
can be solved numerically. As a result, the thicknesses 8" and <5** are 
determined, and then the friction coefficient and Nusselt number are 
obtained: 



C f = (v/US^af) Nu, = (X/ §:%(/, f t ) (1.27) 



where tj(f) and £,(/,/,) are tabulated functions [7]. 

Methods for calculating heat transfer coefficients from nonisother- 
mal surfaces were also developed by Ambrok [11], Shuh [12], Spalding 
[13], and Seban [14]. A greater number of studies, which are surveyed by 
Spalding and Pun [15], are concerned with heat transfer from isothermal 
surfaces. 

Integral methods are extensively used for calculating coefficients 
of friction and heat transfer in turbulent flows. In this case, however, 
experimental data must be used for defining the functions describing 
the velocity and temperature distributions in the layer and for defining 
the relationships contained in the integral equations. Many methods are 
suggested. Most of these pertain to calculating friction coefficients. Much 
less is done on heat transfer from isothermal surfaces, and relatively few 
studies are concerned with heat transfer from nonisothermal surfaces. 
Surveys of such methods are presented, among others, in [1,16-19]. 

One of the simplest methods for calculating friction and heat transfer 
in turbulent flow is based on the power functions for the velocity and 
temperature distributions: 



u 

TI 




T-T, 

T -T 



r y / 7 

y_ 

\ s t j 



(1.28) 



These functions describe rather satisfactorily the distribution of 
these parameters in the turbulent layer, but they do not represent their 
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distribution in the viscous sublayer. Because of this, the skin friction and 
the heat transfer cannot be calculated by differentiating function (1.28), 
as in the case of a laminar layer. In this case, the friction and heat transfer 
coefficients are determined using correlations of numbers of experimen- 
tal data. 

Using the well-known law of the wall u/ u T = 8.74(i/!¿ t /v), the constants 
C = 0.045 and n = -1/4 in the general relation C f = CRe£ may be defined 
as follows [19]. Transforming the first part of Equation (1.28) to variables 
of law of the wall leads to the following expression: 

u/u T = (C/2)- 4 / 7 Re¡ (4,,+1)/7 (i /ü r /v) 1/7 , C f = 0.045 Re^ 1/4 (1.29) 

which should be identical with this law. This occurs only if 4)/ + 1 = 0 
when the Reynolds number exponent becomes zero and (C/2) -477 = 8.74, 
giving Equation (1.29). 

The heat transfer coefficient can be estimated using the fol- 
lowing ratio for the turbulent layer out of the viscous sublayer: 
r/q = A ti (d T /dy)/ y tb (du/dy). Calculating derivatives from Equation (1.28) 
and assuming that turbulent Prandtl number Pr (fl = 1 and that t = x w and 
q = q w in the vicinity of the wall, one gets 

St = 0.0225 Pr _1/4 Re¡ 3/4 (d/<5 ( ) 177 (1.30) 



where Pr _3/4 is added to take into account the Prandtl number effect. 
Using this expression and relation (1.29) for the friction coefficient, one 
obtains from the integral relationship (1.16) two differential equations: 



dd 5/ 4 
dx 



115 1 dU 
28 U dx 



8 5 ' 4 = 0.289 




(1.31) 



de 9/7 t 8 1 d(Ue„8) 0.260 í v] 

dx 9 U9 W 8 dx Pr 374 S 5 ' 4 { U ) 



(1.32) 



Integrating these equations [16] yields 8 and e = S t /8. Then, 
Equation (1.29) and Equation (1.30) can be used for calculating the fric- 
tion and heat transfer coefficients. 

Integral methods based on power law approximation (1.28) as well as 
methods employing polynomials are suitable in those cases when the 
velocity and temperature profiles across the layer do not exhibit singu- 
larities. The law of the wall is valid for a whole boundary layer at small 
Reynolds numbers. At high Reynolds numbers, the velocity distribution 
across the layer is described by a composite function that consists of the 
law of the wall for the inner and the velocity defect law. 



U-u 



= f 



u. 



1 

8 



(1.33) 
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for the outer parts of the turbulent boundary layer. (See details in 
Section 4.1.) Integral methods employing such composite profiles were 
developed by Patankar and Spalding [17], 

In calculating the coefficients of heat transfer from nonisothermal 
walls, extensive use is made of a method based on the conservative 
nature of the laws of heat transfer, written in the form of relationships 
between the local parameters of the boundary layer [20-23]. Theoretical 
analysis and experiments [22,23] show that one can neglect the effect 
of the pressure gradient and thermal boundary condition on the law 
written as 



St = C Re"„ Pr m (1.34) 

This makes it possible to use, in the general case, this equation derived 
for the simplest case of gradientless flow and constant wall temperature. 
Substitution of Equation (1.34) into the second integral relationship (1.16) 
yields an equation that, when integrated, gives 



i 

1 -n 

(1.35) 



The constants C = 0.0128, n = -1/ 4, m = -0.6 are obtained from workup 
of a large volume of experimental data. Once Re á » is defined, St is found 
from Equation (1.34). 

This approach allows determination of the solution of a large number 
of problems in a rather simple manner and with accuracy sufficient for 
practical needs [20,21,24]. However, it is shown [21-23] that the relative 
error in St, resulting from neglecting the effect of boundary conditions 
on the governing relationship, can be estimated by the inequality 



Re,„ = — 

sr q_ 



7 f 

C(l-n)Re L Pr m 9 w 1 ~ n d 

• v 



A St 

17 



< 0.1 



Re C dd w 

StRe L 9 w d(x/L) 



(1.36) 



1.5 Method of Superposition 

The equation of the thermal boundary layer is linear, so one can apply to 
it the principie of superposition. According to this principie, the solution 
for the general case of an arbitrary nonisothermal wall can be obtained as 
a sum of particular Solutions, found for some simple law of the wall tem- 
perature variations. One such simple law is stepwise change in temperature 
downstream of a nonheated zone. In this case, the wall temperature remains 
equal to the free stream temperature to some point x = t, and at this point 
changes suddenly to some constant valué T w . It the solution of this problem 
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is known, then for arbitrary wall temperature distribution, the solution is 
given by the following integral [25,26]: 

X i 

^=- = 6(x,y$, T-T x =\e(x,y^)^d^ + ^d(x r y r QAT wk (1.37) 

W o ’ k=l 

Here the integral pertains to the continuous part, and the summation per- 
tains to breaks in the temperature distribution, and A T wk is a temperature 
jump in a point with coordinate £, k 

Differentiating Equation (1.37) with respect to y and using the influence 
function/(x, q) yields the expression for the heat flux at the arbitrary noniso- 
thermal wall: 



% = K 



X i 

\fM^ + ^f(x,QAT ak 

o ’ fc=i 



h s /h. = f(x,Q (1.38) 



Here the influence function describes the effect of the nonheated zone on the 
coefficient of heat transfer downstream of the sudden temperature change. 
Expression (1.38) is valid for any flow mode, and because successful methods 
for determining h, are available, the main difficulty consists in calculating 
the influence function. 



Example 1.3: Calculating the Influence Function [25,26] 

For the elementary case of laminar gradientless flow and Pr > 1, the 
influence function is calculated by the integral method. The velocity 
and temperature profiles across the layer were approximated by third- 
degree polynomials. In this case, the thickness § and the ratio of e = S t / 8 
are calculated by the same differential equation (Equation 1.16) as when 
using polynomials (Equation 1.18) but with different coefficients: 



5 dS _ 140 v 
dx 13 U ' 



§£ dm =w a 



dx 



U 



(1.39) 



Substituting 8 obtained from the first equation into the second equa- 
tion yields a differential equation with a solution (with 13/14 = 1) that 
gives the influence function in the following form: 

e 3 + 4e 2 x— = — /(£/ x ) = [i_(¿/ x )3/4]-i/3 (1.40) 

dx 14 Pr 

A more general result, pertaining to arbitrary gradient flows but lim- 
ited to high Pr, is due to work by Lighthill [27]. He used the Prandtl- 
Mises form of thermal boundary Equation (1.3) and assumed that the 
velocity in the thermal boundary layer is proportional to y or, according 
to Equation (1.4), to \¡/ 1/2 : u = t w y/ y = (2t k ,i/// y) in . The correctness of this 
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assumption improves at higher Pr, as the thermal layer becomes thinner 
compared with the velocity layer. Substituting the relation for velocity 
into Prandl-Mises Equation (1.3) gives 



dT _ 1 


f 2 t v 


1/2 A 


' ,„2i] 


dx Pr 


l p y 


dy/ 





(1.41) 



Solution of this equation for the heat flux in the case of stepwise tem- 
perature A T w and a varying temperature T w (x) is given by Lighthill [27]: 

^ = 0 . 5384 ^^ fjx) 



Example 1.4: Accuracy of Linear Velocity Distribution [27] 

To estimate the exactness of linear velocity distribution across the ther- 
mal boundary layer, Lighthill calculated the Nusselt number for self- 
similar solution and uniform surface temperature. In this case, Equation 
(1.8) becomes 



X 

y 



Tjz)dz 



A0„ 



(1.42) 




0.807 Pr 1/3 [<p'(0)] 2/3 



3 

2 



1 - 2/3 



(ra + 1) 



(1.43) 



Table 1.2 compares the results obtained by this equation with 
self-similar Solutions for Pr = 0.7. The agreement is best for m = 0, because 
in this case the derivatives of velocity at the wall vanish, and the velocity 
profile is cióse to linear. For m > 0, the error is greater because d 2 u/dy 2 < 0. 
For m< 0, the profile has a point of inflextion, d 2 u/dy 2 is positive for 
small y and negative for large y, giving too small a heat transfer rate. 

Using Equation (1.42) and expression for shear stress, one can cal- 
cúlate the influence function. In particular, for self-similar flows it is 
found [16,26]: 



/(<g/x) = [1 - ( < g/ x )(3/ 4 )( m +l)]-l/3 



(1.44) 



TABLE 1.2 



Accuracy of the Lighthill Formula (Equation 1.42) 



m -0.0904 


-0.0654 


0 


1/9 


1/3 


1 


4 


— » 00 


Equation (1.43) 


0.497 


0.601 


0.648 


0.653 


0.587 


0.398 


-0.921 


Self-similar data 


0.541 


0.585 


0.596 


0.576 


0.495 


0.325 


-0.741 


Error (%) 


-8 


+3 


-9 


-13 


-18.5 


-22.5 


+24 
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Equation (1.43) becomes exact as Pr — > In the second limiting case, 

when Pr — » 0, the thickness of the velocity boundary layer tends to zero, 
so that the longitudinal velocity over the entire thermal layer can be 
taken equal to the free-stream velocity: u = U(x). The equation of the ther- 
mal boundary layer then becomes simpler, and for self-similar flows, the 
influence function can be presented in a form similar to that of Equation 
(1.44): 



... ,c)T d 2 T 

U(X) *c =a W 



/(£/*) = [l-(£/x)’» +1 ]- 1/2 



(1.45) 



The influence functions for turbulent gradientless flow at Pr = 1 are 
similar: 



/(£/*) = [l-(£/x)9/i°]-i/s f £/*) = [ l-M- 0114 (1.46) 



The first expression can be derived theoretically by the integral 
method, using power law distributions (Equation 1.28) for velocities 
and temperatures [26]; the second one is obtained in Reference [21] by 
using the heat transfer law (Equation 1.34) and solving an integral equa- 
tion. Both influence functions in Equation (1.46) yield virtually identical 
results that are in satisfactory agreement with data obtained by several 
experiments [28]. A stronger effect of the unheated zone is described by 
the two other functions [29]: 

/(£/*) = [ 1-(£/x)W«]- 7 /39 f(£/x) = [ l-(^/x)]- 02 (1.47) 

The method of superposition can also be used for solution of the 
inverse problem, when the temperature head distribution d w (x), which 
corresponds to known heat flux distribution, q w (x) should be estimated. 

The first part of Equation (1.38) shows that such an inverse problem for 
continuity temperature head is reduced to the Volterra integral equation, 
with the influence functions for gradientless self-similar flow in the fol- 
lowing form: 



X 

/(£/*) = [ i-(£/*) Ci r C2 ( L48 ) 

Introducing the new variables, z = x Cl ,£ = and F(z) = q m /h t z c 2 , 
one reduces Equation (1.48) to the generalized Abel equation with the 
solution 



Z 

I 



w 

(z-Q c > 



=m 



dT u, = 1 d 

dz (C, -1)(-C,) dz . 






(z-C) 



(1.49) 
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Returning to variables x and £ leads to the equation defining the 
temperature head [26]: 



T-T =- 



C, 

r(C 2 )r(i-c 2 ) 



í- 



Ci(l-C 2 ) 



hm 



di ; 



(1.50) 



1.6 Solutions of the Boundary Layer Equations 
in the Series of Shape Parameters 

The Solutions of boundary layer equations can be represented in the form 
of series in powers of perturbations of homogenous conditions at the outer 
edge of the layer and on the surface of the body. In this case, the first term of 
the series corresponds to the solution of the given problem for homogeneous 
conditions — that is, for constant free-stream velocity U(x) and constant wall 
temperature T w (x). The subsequent terms of the series represent the effect of 
the perturbation of the homogeneous conditions and contain the sequence of 
derivatives of functions U(x), and TJx), (or 0 w {x)) with respect to longitudinal 
coordinates. The series are constructed using dimensionless combinations 
of variables, which, depending on the number of the derivative they contain, 
are called the first, second, and so forth, shapes parameters. 

One uses either sets of shape parameters that depend only on derivatives, 
or those that also correct for the effect of previous history by incorporating 
the thicknesses of the layer. The solution of equations of the velocity bound- 
ary layer in the form of series of shape parameters of the first type 



x k d k U j. _ x k d k 6 w 

Ulhp' tk ~ 0 W dx k 



(1.51) 



is derived in [30]. An analogous solution for the equations of the thermal 
boundary layer is derived in [31]. These Solutions are represented by múltiple 
series containing nonlinear combinations of the shape parameters: 



-r =7 = ¥ 0 (n)+ ^i(»7)/i + V 2 (n)f 2 + ¥u(n)fi + ¥ 3 (l)f 3 + VnWfJi + - ( L52 ) 

\JvxU 



9 = = T 0 (ri) + TfWf, + 7 ^( 7 ?)/,, + Tf ( 77 )/, + T?\rj)f t2 + - (1-53) 

1 W °o 

Here r\ = if-Jll / vx, bottom indices of functions and T(r¡) give the num- 
ber, whereas the upper ones indicate the type of shape parameters (velocity 
or thermal). 
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Substitution of these series into boundary layer equations yields a set of 
ordinary differential equations defining the coefficients V / (b) and T(r¡). The 
friction and heat transfer coefficients (C f /2)Re\ /2 and Nu,Rc\. l/2 are presented 
by the same series in which functions t/(?7) and T(/;) have been replaced by 
the corresponding valúes of derivatives at r] = 0 : t//"(0) and [-T'(O)]. Equations 
defining functions y/(ri) and T(r¡) do not incorpórate the distributions of U(x') 
and T w (x), which vary from one problem to another. As a result, they are 
universal and can be tabulated. 

Analogous results are obtained by using shape parameters of the second 
type, which additionally contain the characteristic thicknesses of the layer: 



fk = 



r £**2 ^ 


i \ T ,,d k U 


^ X **2 ^ 


l v J 


\ u ‘ M 


r 

l v J 



Wd% 
Q,„ dx k 



(1.54) 



In this case, the solution of the boundary layer equations is represented by 
the previous series (Equations 1.52 and 1.53). Only the scales for variables are 
different: U8**/B and 8**/B, where B is a normalizing constant, for t// and i), 
respectively. As in the case when shape parameters of the first type are used, 
the friction and heat transfer coefficients (C / /2)(Re ¿ .„ /B) and Nu,8 ,t /xB, are 
given by series (1.52) and (1.53) in which the functions Yil) and T(r¡) are 
replaced by the corresponding valúes of the derivatives at r¡ = 0 :t//"(0) and 
[-T'(O)]. The results of integration of the corresponding equations defining 
coefficients y/(ri) and T(r¡) and the valúes of derivatives needed for calcula- 
tions are given in the book by Loytsyanskiy [32] and in the paper by Oka [33], 
for the velocity and thermal layers, respectively. 

To be able to calcúlate shape parameters (Equation 1.54), one must first 
find the thicknesses 8" and 8¡* . As the first approximation, this can be done 
by means of the standard integral relationships. Then, <5** and 8" are refined 
using the integral equations based on the share stress and heat transfer flux 
calculated by series (Equations 1.52 and 1.53) [32,33]. The advantage of shape 
parameters of the second kind (Equation 1.54) is that in this case the initial 
thickness of the layer can be arbitrary, whereas when the shape parameters 
of the first type (Equation 1.51) are used, the initial thickness of the layer is 
always zero. This pattern makes it possible to joint different regions of the 
boundary layer using 8** and 8¡*. 

The effectiveness of the shape parameter methods is frequently limited by 
the rafes of convergence of the series. These difficulties can be circumvented 
by reducing the initial equations of the boundary layer to universal form 
[34]. In this case, one can evalúate the universal equations numerically and 
compile tablets that can then be used in specific calculations, without having 
to resort to series. To be able to reduce the boundary layer equations to a uni- 
versal form, one must use the shape parameters as independent variables. 

This approach eliminates difficulties due to series convergence, but other 
problems arise. These are related to numerical integration of differential 
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equations with a large number of independent variables, and also to the use 
of tables with many entries. So far, one cannot expect to obtain Solutions 
of universal equations with not more than three variables [1]. This means 
that for the velocity equation, one can construct either a single-parameter 
solution that is a function of the first shape parameter, or a two-parameter 
solution that is a function of the first two shape parameters. The first of these 
Solutions is equivalent to the exact solution of boundary layer equations with 
a linear velocity distribution U(x) in the free stream, whereas the second is 
equivalent to the exact solution for a parabolic distribution of U(x). Some 
results of calculations are presented by Loytsyanskiy [34]. 
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Approximate Solutions of Conjúgate 
Problems in Convective Heat Transfer 



Beginning in the 1960s, when computers carne into use and it became clear 
that the conjúgate heat transfer formulation was a natural, comprehensive, 
and powerful approach, numerous Solutions of conjúgate heat transfer prob- 
lems were published. Analytical and numerical methods of different levels 
of accuracy were used, however, as not only early works, but frequently con- 
temporary studies as well are based on the approximate analytical methods 
outlined in Chapter 1. 

In this chapter, typical examples of such approximate analytical Solutions 
published during the last 50 to 55 years are presented and discussed. Two 
types of approximate Solutions are considered. The first are Solutions based 
on the approximate energy equation, when, for example, linear velocity 
distribution across the thermal boundary layer or slug flow model is used. 
Another group of problems considered in this chapter ineludes the solu- 
tion for which approximate results could not be refined or were difficult to 
upgrade. Such Solutions are obtained, in particular, by integral methods or 
by using series that slowly converge. 

Exact and more accurate Solutions of conjúgate convective heat transfer 
problems are considered in following chapters. 



2.1 Formulation of a Conjúgate Problem 
of Convective Heat Transfer 



In general, the conjúgate problem of convective heat transfer in the boundary 
layer approximation is governed by a system of four equations: continuity, 
momentum, and energy equations for a fluid and a conduction equation for 
a body. In the case of a two-dimensional heat transient problem for incom- 
pressible fluid, this system is 



du dv 
dx dy 



( 2 . 1 ) 



du du du 1 dp Q d 2 u 

» + “s +B *--?* + ** r - rJ+ v 

dT dT dT d 2 T v f du '' 2 

— + w— — + u — = a- 1- — — 

dt dx dy dy - c p ydy j 



( 2 . 2 ) 

(2.3) 
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dT, i 



( d 2 T ' 3 2 T A 



dx 2 dy 2 



+ C lv 



(2.4) 



The boundary and conjúgate condition are 



y = ° T s= T w K 



y — ^ °° ^ ^ Lí, T — > 


(2.5) 


e-h I 
ro < 

II 


(2.6) 



dy 



v ^ y „= 0 



<H = 0 



The initial conditions are usually specified for each particular problem. 

In the case of torced heat transfer, the second term in the right-hand part 
of Equation (2.2) should be omitted, and the first should be substituted by 
U(dU/dx) = - (1 /p)(dp/dx). This relation follows from Equation (2.2), written 
for the free stream of the boundary layer. In the case of natural convective 
heat transfer, the first term in the right-hand part of Equation (2.2) should be 
omitted. 

There are two cases when the problem for the body can be simplified. If the 
body is thin (A « L) and its thermal resistance is comparable with that of the 
heat transfer agent, (Bi = HA/ A s ~ 1), then the longitudinal conduction can be 
neglected. In such a case, the temperature distribution across a body thick- 
ness is described approximately by the one-dimensional conduction equa- 
tion. That form follows from Equation (2.4) when the derivative with respect 
to the longitudinal coordínate is taken to be zero. If, for example, there are 
no heat sources inside a body, in the case of a steady-state problem, the tem- 
perature distribution across its thickness is cióse to linear, and, henee, there 
is no need to solve the conduction equation for the body. 

Another case of simplification is when the body can be considered as a 
thermally thin object. This means that the body thermal resistance is small, 
so at any instant, the body temperature is practically constant over its thick- 
ness. In such a case, it is reasonable to intégrate the conduction equation (2.4) 
across body thickness to get an averaged form: 



1 dT ap d-T av ^ <7 rol + tj a , 2 (?;,)„» _q 
a dt dx 2 A A A 

s s s 



(2.7) 



2.2 The Case of Linear Velocity Distribution 
across the Thermal Boundary Layer 

There are at least three situations when this assumption is reasonable: when 
a fluid is characterized by a large Prandtl number, in the case of existence 
of an unheated starting length upstream of a heated section, and for the 
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FIGURE 2.1 

The effect of a Prandtl number on the thicknesses of velocity and thermal boundary layers: (a) 
Pr — > 0, b; (b) Pr — > 



entrance región of a tube or a channel with a fully developed velocity pro- 
file. In the first case, the thickness of the thermal boundary layer is small in 
comparison with that of the velocity boundary layer because a large Prandtl 
number indicates that fluid viscosity is large and, henee, the influence of the 
wall goes far into the fluid (Figure 2.1). Consequently, the velocity boundary 
layer becomes thick, and the distribution of velocity across a thin thermal 
boundary layer can be approximately substituted by a tangent. 

In the other case, the thermal boundary layer starts to develop inside the 
velocity boundary layer (Figure 2.2), which grows when it flows along the 
unheated part of a píate. As a result, the thickness of the thermal boundary 
layer is much smaller than the thickness of the velocity boundary layer, and 
the velocity distribution across the thermal boundary layer is cióse to linear. 




FIGURE 2.2 

The effect of an unheated zone on the development of the velocity and thermal boundary 
layers. 
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A similar situation takes place in the entrance región of a tube or a channel 
when the thermal boundary layer starts to develop inside the fully hydrody- 
namically developed flow. 



Example 2.1: An Incompressible Flow with a Large 
Prandtl Number Past Píate with a Heat Source [1,2] 



In the case of a steady-state problem for a píate, the System (2.1) to (2.4) 
has the self-similar solution with ¡5 = m = 0 (Section 1.1). In this case, 
according to Equation (1.6) and Table 1.1, the li near v elocity profile can be 
presented as (p(r¡) = cr¡ = 0.664 rj, where cp = i y/\lvUx and r¡ = (y/2)y¡U/ vx. 

The energy equation for fluid is used in the Prandtl-Mises form 
Equation (1.3), and the conjúgate condition (2.6) is presented as follows: 



dT d f dT) 

^- a ^{ u ^r 



dT 

ay/ 



y /-0 



= -AS(x) 



( 2 . 8 ) 



where A = (A S /A) and S{x) = -{dT /dy} 0 . Using these two relations and 
a new variable z = ( 2/3)y/ 3/i , one obtains a differential equation with 
boundary condition 



4 i/ 4 9T_a 2 T j^ar z 1 ' 3 ( ar i 

b dx dz 2 3z dz 2 4/3 a x 1/4 1 dz I n 

\ y z-0 



= -AS(x) 



b _ a^Jc/2U 3 ' 4 
v 1/4 



(2.9) 



where the fluid temperature far away from the píate is taken to be zero, 
and c is constant in the velocity distribution. Equation (2.9) can be trans- 
formed to an integral equation [1,2]: 



T(x) = 



3 1/3 aA f S(y) 
2r(2/3)Í7 2 / 3 J (x 3 ' 4 -y 3 ' 4 ) 2 ' 3 Uy 



( 2 . 10 ) 



The solution of conduction equation (2.4) for the píate is obtained by a 
standard method. For example, in the case of the heat source (Equation 
2.11), one obtains 



1 x>0 
0 x<0 



T(x) = Je(Z-x) 



q(x>y)= c i A l ~ x ) £ ( x )= 

^ | G(x,y)S(i/)dy 



( 2 . 11 ) 
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The System of two integral equations (Equations 2.10 and 2.12) deter- 
mines two unknown functions T(x) and S(x). The asymptotic solution of 
this System for x/A » 1 is obtained by Perelman [1,2], The first terms of 
the series defining the surface temperature and the Nusselt number are 



T(x) = 2V3 q X 

’ A. 3 2/3 c 1/3 S 



1 1 ~ A 


1 

- 

IN 


2 5/3 3 2/3 c l/3 ^ x 


5 . 2 i/3 3 V3 c m s [ x ) J 



(2.13) 
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Re 1 / 2 



2 1/3 ^ 
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r(2/3) 



(2.14) 



Re V2 p r l/3 



81 = 



Re 1 / 2 Pri/ 3 



These results are applicable for (a/ Re 1 / 2 Pr 1/3 ) « x < l , Pr > 1, Re^ < 3 ■ 10 5 . 

The same method of solution is used in two other early studies. Kumar 
and Bartman [3] considered the same problem for compressible flow 
over a radiating píate. Kumar [4] solved the problem for a semi-infinite 
porous block with injection streamlined by an incompressible fluid. 



Example 2.2: Transient Heat Transfer in the Thermal Entrance 
Región of a Parallel Píate Channel with the Fully Developed Flow [5] 

First, a simple problem for a channel of height 2 R is considered using 
mass average velocity u m (slug flow). A flow and channel are initially at 
constant temperature T 0 when suddenly the temperature of the channel 
walls becomes an arbitrary function of position and time, and the fluid 
inlet temperature begins varying arbitrarily with time. The governing 
equation and boundary conditions are 



30 

dt 



+ ll 



m 



dO 

dx 




(2.15) 



f = 0, x>0, 0 <y<R, 0 = 0, x = 0, t>0, 0 <y<R, O = 0 e (t) 
f>0, x>0, y = R, 30/3y = O, y = 0, 8 = 6 w {x,t) 

Here, 0 = T - T 0 and 0, are the inlet temperature excess. 

The next simplification is the use of the quasi-steady approach when it 
is assumed that the steady-state relations are valid at each instant of time 
(Section 6.6). Physically, this means that the heat transfer process in the 
fluid is slow in comparison with that in the walls. The previous System 
of equations becomes 



{T = t-x/uJ 

30 3 2 0 

ii m — = a—— 



(2.16) 
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x = 0, O <y<R, 6 = d,{ r), x>0, y = R, d9/dy = 0, 

y = o, e=e w (x,T+x/uj 

This problem is solved by the Laplace transform to obtain the expres- 
sion for surface heat flux using a new variable a = t - x/ u m + £/ u m : 






Y (-1)" exp{-[R 2 n 2 / a(t - o - )]} 



y] Ka(t-o)/R 2 



T,o)do for t<x/u m 
da 



M-iJ* 



¿ (-1)" exp|-[R 2 « m w 2 /a(x - 1)]} 



j7ia(x-Q/u m R 2 



fe, <2 ' 18) 

o? 



It is asserted that these equations will hold approximately in nonslug 
flows if the kernels in the integrands come from the appropriate nonslug 
flow Solutions. Thus, if a solution for the surface heat flux in steady, non- 
slug flow over an isothermal surface is 

q Jx) = G(x)e w0 , (2.19) 



where 6 m0 is a constant wall temperature excess, then Equations (2.17) 
and (2.18) can be written for t < x/ u m and t> x/ u m , respectively, as 

í* r)9 í* df) 

<?„, = J G[u m (t-a)]-^(T,a)da q w = J G(x-q) (^T)d^ (2.20) 

o o s 

Expression (2.20) take into approximate account both thermal history 
and finite thermal capacity of the fluid, while other authors who used 
the quasi-steady approach focused only on the effect of thermal capacity 
[6] or on the effect of thermal history [7,8]. 

Sucec [5] used expression (2.20) to solve two transient heat transfer 
problems. He considered a fully developed flow in a parallel píate duct 
with outside insulated walls. In the first case, a transient is initiated 
by step change in the fluid inlet temperature; in the second case, the 
unsteadiness is caused by a fluid inlet temperature that varies sinusoi- 
dally with time. For a thermally thin píate, the energy balance is used in 
the form that follows from Equation (2.7): 

=Q 

Pw C piA 



( 2 . 21 ) 
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The linear velocity profile is chosen to approximate the actual velocity 
distribution in the entrance región of a duct with fully developed flow. 
Using the mass average velocity, u m , and Lighthill's formula (1.42) for the 
steady-state heat flux, one obtains function G{x) in Equation (2.19), where 
L = ca/ u m R 2 is the nondimensional distance along a duct: 

u x (y) = 3u m (y/R) G(x) = 0.5384(3) 1/3 (A/R)L _1/3 (2.22) 

If T¡ is an initial temperature, 0 w = T w - T, and henee, = 0 at 
t = 0, then it is seen by inserting the first part of Equation (2.20) into 
Equation (2.21) that 0 = 0 for the domain t<x/u m . Inserting Equation 
(2.22) into the second part of Equation (2.20), and then Equation (2.21) 
gives the equation for surface temperature distribution, after noting 
that dO w /dt = dQ m /dx; 



d0 w 0,538(3) 1/3 A f 1 90, , 

p u¡ c pw RA J (£-§)V3 ^ 



(|, x )rf| = 0 for r > 0 



(2.23) 



Applying the Laplace transform to solve this equation subjected to 
side conditions 



x = 0, L > 0, Q w = 0, L = 0, x > 0, 0 w = 6 e 
gives the wall temperature and heat flux for x > x/u ln : 



(2.24) 
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(2.25) 



1 z 



AKZ) ~ k\ 3 Ki/3 { 3 



— 7 3/2 



/ = 



0.538 1/3 r(2/3)Ar 
p c A RÜ /3 
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where Ai(z) and K 1/3 are, respectively, the Airy function and the modi- 
fied Bessel function of second kind. For domain t < x/ u m , the surface heat 
flux and both temperatures are zero. 

The other problem is solved by the same procedure. Equation (2.23) 
is integrated subjected to the same side conditions with 9 e = AT 0 sinft)f, 
where AT 0 is the amplitude of the periodical inlet temperature. The 
results are given in [5]. 

To ascertain the accuracy and limitation of the suggested method, the 
obtained results are compared to the numerical Solutions. It is found that 
an improved quasi-steady approach, which attempts to take into account 
both thermal history and thermal capacity effeets, prediets the finite- 
difference results reasonable well and does significantly better than the 
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standard quasi-steady method. Yet the improved quasi-steady approach 
is more difficult to apply than the standard quasi-steady solution. 

A linear velocity distribution across the thermal boundary layer is also 
used for studying heat transfer between two fluid flows separated by a 
píate without thermal resistance [9]. 



2.3 The Case of Uniform Velocity Distribution across 
the Thermal Boundary Layer (Slug Flow) 

The other limiting case corresponds to a fluid with small Prandtl number. In 
this case, the thermal boundary layer is thick in comparison with the veloc- 
ity boundary layer because a small Prandtl number indicates that fluid vis- 
cosity is small. Henee, the influence of the wall is confined to a región of fluid 
cióse to the wall, and the velocity across the flow is approximately equal to 
the free stream velocity (Figure 2.1). For some fluids, like mercury and other 
liquid metáis, such an assumption is cióse to reality. 

Another reason to consider the velocity components to be constant across 
the thermal boundary layer is to simplify the energy equation. Because of 
the relative simplicity it affords, this assumption is used by a number of 
investigators. Perelman [10] and Rizk et al. [11] studied conjúgate heat trans- 
fer to laminar flow from heated blocks, Luikov [12] investigated the case of 
a píate with a heat source, and Viskanta and Abrams [13] used the slug flow 
assumption to study heat transfer between two concurrent and countercur- 
rent streams separated by a thin píate. 

Several researchers used the slug flow approximation to study transient 
coupled heat transfer when the unsteadiness is caused by different sources: 
by flow with time-varying temperature, by step temperature change of free 
stream at the leading edge of a píate, and by unsteady sources in the wall. 
Sparrow and DeFarias [14], Siegel and Perlmutter [15], and Namatame [16] 
considered coupled heat transfer in the channel flows, and Sucec [17-19] and 
Chambre [20] solved conjúgate problems of heat transfer to flows over piafes. 
Travelho and Santos [21] studied the effect of wall heat capacity in a circular 
duct for the case of periodically varying inlet temperature of the laminar 
flow. 



Example 2.3: Transient Heat Transfer from a Píate of Appreciable 
Thermal Capacity with Heat Generation of q 0 to Turbulent Flow [6] 

The governing equation is the energy equation (2.3) with additional 
eddy diffusivity e H that is estimated using the expression of frictional 
stress: 



dT dT 
— + «— 
dt dx 
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dy 



(a + e H ) 



dT 
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0.0228 Re¡ 1/4 



(2.26) 
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and 1/7 power law (Equation 1.28) of velocity distribution. Assuming 
that turbulent Prandtl number equals unity and using an approach simi- 
lar to that described in Section 1.4, total momentum and heat diffusivity 
are obtained as 



a + e H = 0.1 791Í 27/35 v 8/35 x _3/35 (2.27) 

Because the píate is assumed to have appreciable thermal capacity but 
negligible thermal resistance, the initial, boundary, and conjúgate condi- 
tions are as follows: 



t = 0, x = 0, y — > T = 0, lim 

y^>0 



-pc(a + e H ) 



9T 

3j/_ 



1 i c)T i 

= %~2 P- c - A t>0 ' x>0 

V /y=0 



(2.28) 



For early times when the thickness of the thermal layer is thinner than 
that of the laminar sublayer, the eddy diffusivity can be neglected in 
comparison with the molecular diffusivity. In this case, the transient 
energy equation (2.26) reduces to the laminar slug flow equation that has 
a solution for the case of a step heat source that has two forms: for small 
and at the far downstream position where ü < X, and for large and at the 
leading (upstream) position where d > .X, respectively: 

(2.29) 



(2.30) 
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Here, ú = t/B 2 , X = x/B 2 U, T RS = c¡ 0 4aB/X is reference temperature; 
B = cVce / A is the thermal capacity parameter; and c = (l/2)p n ,c n ,A is píate 
thermal capacity. These results are valid for times less than 1 ms. For 
times greater than 1 ms, the thermal boundary layer penetrated into 
the región represented by the 1/7 power law (1.28). The solution of tran- 
sient turbulent energy equation (2.26) together with Equations (2.27) and 
(2.28) is found by two steps considering several simpler problems. The 
first is the case of zero thermal capacity. Two equations are solved for 
this case: 
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dt 
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(2.31) 



The first equation corresponds to small but greater than 1 ms times, 
for which convective terms have little effect. This equation is solved 
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with initial condition t = 0,T = 0. The second equation (2.31) corresponds 
to large valúes of time when transient effects should die away, and the 
heat transfer characteristics should approach those of a steady-state con- 
dition. Consequently, the initial condition for this equation should be 
x = 0,T = 0. Introducing a new variable x/U, one can see that in the case 
of the slug flow, the second equation (2.31) is the same as the first. Thus, 
both equations have similar Solutions: 



pcUTJt,x) 

% 

p c UTJt,x) 

% 



= 30.3 Re 0 - 2 



tu ^ 
' T * = lc ~ 1 



(2.32) 



= 30.3 Re° x 2 , t >1 



Averaging these equations along the píate yields the mean wall 
temperature: 



Ut) = 



q ° L (28. 19r°- 125 - 2.94-r 1 - 2 ) 
A Re»- 8 




TJt) = 25.25 



A Re»- 8 



T>1 



(2.33) 



In the case of appreciable thermal capacity, part of the heat gener- 
ated is stored within the wall, and the rest is transferred in the fluid. 
To determine these parts of heat, one should take into account that the 
corresponding heat transfer equations are coupled. Integrating along the 
píate an energy balance at the wall and assuming that heat at the wall is 
proportional to heat generated, q = q 0 f(t), resulting in the expression for 
function /(t): 



q — q 0 c 
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m = i 



i gw 

X L dz 



(2.34) 



Rearranging Equation (2.33) using Equation (2.34) leads to first-order 
differential equations [6]: 



d< J m / T Rs) + Re t 8 x l C jt ) = X L r<l, F(t) = 28.19f° 125 - 2.94-r 1 - 2 (2.35) 

dz F(z) 

d(TJT Rs) + RerX L C ( 2 . 36 ) 

dz 25.25 m RS L 

where C = B\fa /L = ca/X is the dimensionless thermal capacity param- 
eter. Equation (2.35) is solved for the mean wall temperature subject to 
initial condition T m (0)/T RS = 0. This solution at z = 1 is then used as the 
initial condition for Equation (2.36). The final results for the mean tem- 
peratures are 
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(2.38) 



The mean heat transfer coefficient for the píate of appreciable capac- 
ity is defined by relation, which after using Equations (2.37) and (2.38) 
yields 



_ <? 0 -cjdTJdt) 

m T ' 

1 m 



A Re» 8 
LF(t) ' 



A Re 0 - 8 

t<1 h m = 0.0396^-^- 



r > 1 



(2.39) 



Analytical results for the wall temperatures are compared with 
authors' experimental data for two valúes of the dimensionless wall 
capacities (¿ = 0.68 • 10 -3 an d 2.39 ■ 10“ 3 ) and for Reynolds number range 
of 5.3 • 10 5 < Re L < 2 ■ ! O 6 . Agreement of turbulent flow solution (Equations 
2.37 and 2.38) with experimental data for all these parameters is observed. 
The slug flow solution (Equations 2.29 and 2.30) overestimates the wall 
temperatures. The experimental results show that the transient mean heat 
transfer coefficient for high wall thermal capacity and high flow velocities 
decreases gradually until it passes through a mínimum and then starts 
rising, finally reaching its steady-state valué. However, the valúes of mean 
heat transfer h m obtained from the solution do not exhibit any mínimum, 
so that the steady-state valué of h m is reached abruptly at r = 1. 

Kawamura [22] used a similar approach to take into account the ther- 
mal capacity in the case of turbulent flow in an annulus. 



Example 2.4: Unsteady Heat Transfer between a Fluid 
with Time-Varying Temperature and a Píate [17,18] 

A thermally thin píate with lower surface temperature is exposed to a 
coolant at temperature T with heat transfer coefficient h c , which is stream- 
lined over its top by steady laminar flow with slug velocity distribution. 
Initially, the píate and fluid are both at the coolant temperature T, when 
suddenly the fluid temperature at x = 0 is changed to T 0 and subsequently 
held constant. Defining 9 = T - T 0 , one gets the energy equation (2.3) and 
boundary and conjúgate conditions (2.5) and (2.6) for the case of transient 
heat transfer and slug velocity in the boundary layer as follows: 
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Conjúgate Problems in Convective Heat Transfer 



t = O, x > O, y > O, 9 = 9 c = T c -T 0 , x = 0, t > O, y > O, 9 = 0, 
y — > 09 t > O, 0 is /mite 



(2.41) 



y = 0, t>0x>0 c^-+-j-(9-9 c ) = ^- (2.42) 

at Á. ai/ 

where c = p w c w A/A. The problem is solved using a Laplace transform 
with respect to both independent variables. 



9 = J 9exp(-pt)dt 9 = J 9exp(-sx)dx 

o o 

The System of Equations (2.40) through (2.42) becomes 
d 2 9 Us + p 0 



(2.43) 



dy 2 



y — > 9 is finite y = 0, 



cW 

dy 



9 = — L 
sa 



-cp 



(2.44) 



0-0 



c h c 
s spÁ 



The solution of this ordinary linear equation and the inverse trans- 
form yields [18] 



0/ 9 C = 1 - o(t - l){erf(Y) + exp(2 rjY + r 7 2 ){erf[e(r-l) + r¡ + Y] - erf(?j + Y)]] 

Í0 t< 1 (2-45) 

T>1 



o(r-l) = 



Here, z = Ut/x, Y = (y / 2)^ju / ax , r¡ = ( h c / X)^jax/U , e = yJx/aU /2c is 
the coupling parameter, and cr(r-l) is the unit step function. Putting 
Y = 0 in Equation (2.45) gives the píate temperature and then the dimen- 
sionless heat flux scaled by that for the isothermal surface: 

0 n ,/0 c = 1 - c(r - l)|exp(í? 2 ){erf[e(f - 1 ) + 77 ]- erf(^)}} (2.46) 



_ 

A.0 c Vlz7 nax 

\ 0 T < 1 

|exp{-[e 2 (T - 1) 2 + 2r¡e(z - 1)]} + 4ñr]exp(ri 1 ){erf[e(-r - 1) + r¡] - etí(ri)} z> 1 



(2.47) 




